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Abstract 

A method is presented to investigate diffraction of an electromagnetic plane 
wave by an infinitely thin infinitely conducting circular cylinder with longi- 
tudinal slots. It is based on the use of the combined boundary conditions 
method that consists on expressing the continuity of the tangential compo- 
nents of both the electric and the magnetic fields in a single equation. This 
method proves to be very efficient for this kind of problems and leads to fast 
numerical codes. 

I. INTRODUCTION 

The problem of the penetration of electromagnetic waves in a conducting circular cavity 
through a narrow axial aperture has been treated by several authors. Several methods have 
been used to achieve the determination of the field inside the cavity. Beren [Q used the 
Aperture Field Integral Equation , the Electric Field Integral Equation and H-field Integral 
Equation to determine the field around an axially slotted cylinder, while Johnson and Zi- 
olkowski gave a generalized dual series solution for this problem. Mautz and Harrington 
treated the field penetration inside a conducting circular cylinder through a narrow slot in 
both TE || and TM 0] polarizations. More recently Shumpert and Butler ||, |J proposed 
three methods to study the penetration in conducting cylinders. In this article, we propose 
a method to calculate the field inside and around a slotted circular cavity with longitudinal 
slots. It is based on the combined boundary conditions method introduced first by Montiel 
and Neviere Section II is dedicated to the description of the theory. In section III we 
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give some details about the numerical scheme and then compare our results with previous 
work. 



II. THEORY 



The structure under study is depicted in Fig 1. The space is divided into two regions, region 1 
(exterior region : r > R) and region 2 (interior : r < R) that are assumed to be dielectric and 
homogeneous with relative dielectric permittivities E\ and 62, respectively. On the interface 
between these two media are deposited a finite number of infinitely conducting, infinitely 
thin circular strips that are invariant along the z direction. The device is illuminated by 
a TMz (electric field parallel to the z axis) or TEz (magnetic field parallel to the z axis) 
monochromatic electromagnetic wave under incidence 9 with vacuum wavelength A. 
Throughout this paper we assume an exp (—iut) time dependence. The z component of the 
electric or the magnetic field will be denoted by u (9, r) . We denote by fli the union of the 
strips and by Q2 its complementary in [0, 2tt] . 
In the exterior region we express the total field as : 

ui (6,r) anJ n (kir)exp(inO)+ b n H^ (for) exp(m0) (1) 

Likewise in the interior region the total field may be expressed as : 

U2 (9, r) =y~] c n J n (k 2 r) exp(m6>) (2) 

where a n , b n and c n are the amplitudes of the incident, the diffracted and the transmitted 

waves respectively. We denote J n and Hn^ the Bessel and the Hankel functions of the first 
_ 2tt 

kind. k v = ko^Je^ = —r--s/£p~, with p — 1, 2 and Z denotes the usual set of relative integers. 

A 

Amplitudes a n being known, the problem is to determine amplitudes b n and c n from which 
the total field can be calculated everywhere. For that purpose one must write the bound- 
ary conditions at the interface between both dielectric media. This is done in the next 
subsections by distinguishing the TMz and the TEz cases of polarization. 



A. TMz polarization 

The boundary conditions applied to the tangential components of the electromagnetic field 
at the interface defined by r = R lead to : 



mi (6, R) = u 2 {9, R),We [0, 2tt] (3a) 

( dr ) (e,R) ( dr ) (e,R) ,WeQ2 (3b) 
With the supplementary condition that the electric field must vanish on the strips: 

mi (9, R) = u 2 (9, R) = 0, V9 G fii (4) 
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Following Montiel and Neviere |7|||, equations (|3"15| ) and @ can be combined in a single 
equation that holds for every 8 in [0, 2n] : 



f du 2 \ 
\~dVJ 



( du\ 



(e,R) 



\ dr 



(9,R) 



0, V0 G [0, 2tt] 



(5) 



(1-X(0))u 2 (9,R)+g X (e) 
where x (#) is the characteristic function of set 

and $ is some numerical parameter introduced for dimensional and numerical purposes. 



(6) 



1 if x eVL 2 
elsewhere 



Remark that the set of Eqs. (|3a|), (|3b|) and (f|) is equivalent to the set of Eqs. (|3a|) and 
Since x ($) is 27r-periodic it can be expanded in Fourier series: 



Reporting equations ([]]) and (0) into equation fl3a|) and projecting on the (exp(m6*)) ngZ basis 
gives : 



a n J n {hR) + b n H^{hR) = c n J n (k 2 R), Vn G Z 



(7) 



then reporting equations ([[]), (@) an d ® i n t° equation (|5|) and projecting on the 
(exp(m6*)) ngZ basis leads to : 



c n J n (k 2 R)- Xn~pC p Jp(k 2 R) + 



9 E Xn-p k 2 {c p J p (k 2 R)) - h (a v J'{k 1 R) + bpH^XhR) 
pel, 1 v 



0, 



Vn G Z 
(8) 



where the primes denote derivation with respect to r. From Eq([7|) one can extract c n : 



J n (k 2 Rf n+ J n (k 2 R) 



Vn G Z 



(9) 



and report its expression into Eq(|]) to obtain the following linear system linking the ampli- 
tudes a„ and b„ : 



a n J n (kiR)+ E 

p&L 

-b n H<P(hR)+ E Xn-pb P 



-J P (kjr) - gk 2 ^^-J' p (k 2 R) - ghd'ihR) 



H^ihR) - gk 2 



J P {k 2 R) 

4 1] (hR) v 



J P (k 2 R) 



J'Jk 2 R) + ghH^XhR) 



(10) 



The solution of the linear system ( |T0| ) gives the unknown amplitudes b n and then Eq.(Q) 
gives the amplitudes c„. Thus the field can be computed everywhere in space using Eqs.(|]) 
and ®). 
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B. TEz polarization 



For this case of polarization, the continuity of the tangential components of the electromag- 
netic field at the interface defined by r = R leads to : 



dui 



Si V dr 



dUn 



rhr =- X , We [0,2* 



£2 V ^ r 



mi (9, R) = u 2 (9, R), We Q 2 
With the supplementary condition that the electric field must vanish on the strips: 



du\ 



£1 V dr 



\ ( du 



- -rH =- -H 1 = o,V0eOi 



(e,R) 



£2 V ^ 



Here again we can replace equations ( |llb| ) and (|12| ) by : 



;i-xW) 



du 2 
£2 V rf?" 



+ (0) [«2 (9, R) - Ul (9, R)} = 0, V9e [0, 2tt] 



(11a) 
(lib) 

(12) 
(13) 



(0,R) 



Reporting equations (fj) and (fj) into Eq. (|lla|) and projecting on the (exp(m#)) neZ basis 
gives : 



a n J' n (hR) + b n H^'(hR) = ^ £ -± Cn J' n {k2R)i G Z 

fcl£ 2 



(14) 



Remark that the set of Eqs (|lla|) , ( |llb|) and (fT2f) are equivalent to the set of Eqs (|llb| ) and 
(|T3|). Reporting equations (0),(0) and (0) into Eq.(|I~3"D and projecting on the (exp(m#)) n6Z 
basis leads to : 

—c n J' n (k 2 R) Xn- P CpJp(k,2R)+ 

£ 2 £2 P ez 



9 E Xn-p c p J p (k 2 R) - (a p J p (kiR) + b p H^\k 1 R) 



VnGZ 



(15) 



From Eq.(14) one can extract c n : 

ki £2 ( J'rXhR) , H^XhR) 



k 2 £i \J^k 2 R) 



a n + 



Vn G Z 



(16) 



and report its expression into Eq.(|15|) to obtain the following linear system linking the 
amplitudes a n and b n : 



ki 

a n — J' n (kiR)+ Xn-pdp 
£ i p& 

-b n —H ( n ) '{kiR)+ J2 Xn-pbp 
£ i pel 



-^J p (hR) + 9^^§MhR) - gUhR) 

£l K-2 £l Jp\ K 2ft) 

k ^'(k lR ) - 9^^rH^( kl R) + g&faR) 
£1 h e x J' p {k 2 R) 



(17) 



The Solution of the linear system ( |TTD gives the unknown amplitudes b n and then Eq.(|i~6D 
gives the amplitudes c n . Thus the field can be computed everywhere in space using Eqs.(|]) 
and ©. 
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III. NUMERICAL RESULTS 



The infinite linear systems ( [TDD anc ^ (HIP are truncated to a finite size by retaining only 
(2N + 1) coefficients and solved to obtain a representation of the field at truncation order 
N. The convergence of the results has been checked by increasing integer N and using the 
usual criteria of energy balance ( optical theorem ) and reciprocity. We have also verified 
that the boundary conditions are fulfilled, for instance the nullity of the tangential electric 
field on the strips. In all the calculations carried in this paper we set g = — 10~ 3 . However, 
as mentioned in ||, numerical experiments show that only the sign of g is of importance: the 
numerical scheme is more stable with a negative value of g. All the computations reported 
have been obtained on a Personal Computer (200 MHz processor with 32 Mo of RAM), only 
a few seconds are necessary to perform each result shown here. 

In the following we provide some numerical examples and compare our results with those 
obtained in previous works ||, @ and ||. 

In our first example, we consider a circular cavity with a single longitudinal narrow slot 
(see Fig. 2) with cf> = 5° and we compute the interior field on x axis. Figures 3 (a) and (b) 
show the magnitude of the normalized electric field in both the TM Z and the TE 2 cases of 
polarization. It can be seen that our results are in excellent agreement with those published 
recently by Shumpert and Butler ||, ||, see for instance figure 6 in this last reference. In the 
second example, we consider a circular cavity with an aperture such that = 5°. In figures 
4 (a) and (b) are plotted the normalized electric field amplitude at the center of the cylinder 
for various values of the parameter k\R for both the TM Z and the TE 2 cases of polarization. 
These curves agree with those obtained by Mautz and Harrington (see references Jl0| and 
|TT| ) . It is worth noticing that the resonances in these plots correspond to the modes of the 
cavity. 

Finally we give the map of the electric field around and inside the slotted cylinder when 
excited by a plane wave such that k\R corresponds to a mode of the closed cylinder. We 
can see in figures 5 (a) and (b) that the modes TM i(fcii? = 2.404) and TMn(kiR = 3.832) 
are excited inside the structure. 



IV. CONCLUSION 

We have developed a very efficient and fast method adapted to study diffraction of an elec- 
tromagnetic wave by a finite number of infinitely thin, infinitely conducting strips deposited 
on a dielectric cylinder. It is based on the combined boundary conditions method. The 
method is very low CPU-time consuming. The numerical examples that have been given to 
illustrate the method are not restrictive. One can use as an incident radiation a beam of any 
shape. It suffices to calculate its corresponding incident amplitudes a n . It is also possible 
to study the radiation pattern of a source located at the center of the cylinder by making 
slight changes in the equations. 
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Figure captions 



Figure 1 : Geometry of the problem : a TEz or a TMz polarized plane wave illuminates 
the slotted cylinder 

Figure 2 : Electromagnetic penetration into a circular cavity through a narrow slot. 
Figure 3 : (a) Magnitude of normalized electric field on the x axis of a slotted circular 
cylinder excited by TM Z plane wave {k\R = 0.7, 9 = 180°, <fi = 5°} 
(b) Magnitude of normalized electric field on the x axis of slotted circular 
cylinder excited by TE Z plane wave {k\R = 0.7, 9q = 0°, <p = 5°} 
Figure 4 :(a) Normalized electric field amplitude at the center of the cylinder for various 
k±R with (0„ = O°,0 = 5°) 
(b) Normalized electric field amplitude at the center of the cylinder for various 
kiR with (# = O°,0 = 5°) 
Figure 5 : Map of the electric field for values of k\R corresponding to the modes : 
(a) TMoiand (b) TM n of the cylinder. 
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